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Abstract. We show that a bacteria and bacteriophage system 
with either a perfectly nested or a one-to-one infection network 
is permanent, a.k.a uniformly persistent, provided that bacteria 
that are superior competitors for nutrient devote the least to de¬ 
fence against infection and the virus that are the most efficient at 
infecting host have the smallest host range. By ensuring that the 
density-dependent reduction in bacterial growth rates are indepen¬ 
dent of bacterial strain, we are able to arrive at the permanence 
conclusion sought by Jover et al [Tj. The same permanence results 
hold for the one-to-one infection network considered by Thingstad 
[7] but without virus efficiency ordering. Additionally we show the 
global stability for the nested infection network, and the global 
dynamics for the one-to-one network. 


1. Introduction 

Jover, Cortez, and Weitz [I] observe that some bipartite infection 
networks in bacteria and virus communities tend to have a nested struc¬ 
ture, characterized by a hierarchy among both host and virus strains, 
which determines which virus may infect which host. They argue that 
trade-offs between competitive ability of the bacteria hosts and de¬ 
fence against infection and, on the part of virus, between virulence and 
transmissibility versus host range can sustain a nested infection net¬ 
work (NIN). Specifically, they find that: “bacterial growth rate should 
decrease with increasing defence against infection” and “the efficiency 
of viral infection should decrease with host range”. Their findings 
are based on the analysis of a Lotka-Volterra model incorporating the 
above-mentioned trade-offs which strongly suggests that the perfectly 
nested community structure of n-host bacteria and n-virus is perma¬ 
nent, or uniformly persistent M- 
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Inspired by their work, in [6] we replace the Lotka-Volterra model 
by a chemostat based model in which bacteria compete for a growth- 
limiting nutrient. In a chemostat model, each bacterial strain is en¬ 
dowed with a break-even concentration, R*, of nutrient below which 
it cannot grow such that, in the absence of virus, only the strain with 
smallest R* survives. Thus, within a community of bacteria competing 
for a single limiting nutrient, the competitiveness of the various strains 
are naturally ordered by their R* values. In [6], we show that nested 
infection networks of arbitrary size are permanent provided that R* 
values increase with increasing defence against infection and that the 
efficiency of viral infection should decrease with host range. We also 
show how a bacteria-virus community with NIN of arbitrary size can 
be assembled by the successive addition of one new species at a time, 
answering the question of “How do NIN come to be?”. 

We show that the Lotka-Volterra based model of Jover et al [lj can 
be modified in such a way that the permanence conclusions which they 
sought can be attained. The key is to ensure that density-dependent 
reduction in bacterial growth rates be independent of bacterial strain. 
Following pp, we assume that virus strain V{ is characterized by its 
adsorption rate </>* and its burst size /3;, both of which are assumed to 
be independent of which host strain it infects, and its specific death rate 
The density of bacteria strain i is denoted by B i: and its specific 
growth rate is r*. The “mean field”, density-dependent depression of 
growth due to inter and intra-specific competition term JV djBj is 
common to all strains. The equations of our model are the following. 



n 



where matrix M captures the infection network structure: 



1, Vj infects H L 
0, Vj does not infect Hi 


In the system considered in [1], the bacterial host dynamics in the 
absence of virus is modeled as B[ = 7yl?,;(l— K _1 JV Bj); a consequence 
of this is the simplex of equilibria JV Bj = K. We avoid this. 
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Motivated by the work of Jover et al |Tj and the work of Thingstad 
[7], we consider two special network structures: nested infection net¬ 
works (NIN) with upper triangular matrix M, and one-to-one infection 
networks (OIN) with M — I, the identity matrix. 

The scaling of variables 


Pi = <f>iVi, Hi = Bi, ei = 


Pi Oi 


rii 

exposes a virus infection efficiency parameter e* for each virus, 
after, we consider the resulting scaled system: 


Here- 


( 1 . 2 ) 


ft,' = ft 


3 = 1 


P'i = etriiPi ( ^2 M jiHj 

o'=i 


ft E 

3 = 1 

A 

- ] , 1 < i < n. 


System (II.2p defines a dissipative dynamical system on the nonneg¬ 
ative orthant of R 2n 


Proposition 1.1. Solutions of (ll.2[) with nonnegative (positive) initial 
data are well-defined for allt > 0 and remain nonnegative (positive). In 
addition, the system has a compact global attractor. Indeed, if F(t) = 


n 


n 


Y Hfit) + Y 


pm 


then 


Fit) < S + (F(0) - S)e~ m < max{F(0),^}, 

and 

n 

linisup F(t) < ^(1 + 
t-y oo “ W ai 

i=i 

n ra 

where K = max {//,;(()), — W = ruin In,;} and Q = YlW + rf)K. 

i=l w ai i=l 


Proof. Existence and positivity of solutions follow from the form of 
the right hand side. Therefore, Hfit) < Hfit)(ri — aiHfit )). Hence 
Hfit ) < K and lira sup^^ Hfit) < ri/ai. 


dF_ 

dt 


Er.ft - (E ft)(E%ft) - E A E'-.ft 

1=1 1=1 j — 1 2=1 * 1=1 


= ^(W + r l )H l -WF. 
1=1 


W Y. 


i=1 


emi 
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The estimate on F(t) follows by bounding the first summation by Q 
and integrating; the estimate on the limit superior follows from the 
estimate of the limit superior of the Hi above and by integration. □ 

2. Nested Infection Networks 
If M is upper triangular, then our system becomes: 


( 2 . 1 ) 


H[ = hA n- £ p. 


3 = 1 


j>i 


P' = e, 


i n,P, H i ~ j ■ 1 


< i < n. 


Let E* be the equilibrium of the system where each component is 
positive. From the second equation of (12. ip . 

j>i 


ei 


e ? _i 


These are all positive if 

(2.2) ei > e 2 > e 3 > ■ ■ ■ > e n 

n 

Let Fi(H ) = r t — a,jH,j , then at E* from the first equation of (12.11) . 
Fi{H*) = E PI Let 1 

j>i 

(2.3) = A + (Oi _ OJ) + (Oi _ Oi) + -h (— _ -±!_) 

ei e 2 ei e 3 e 2 e n e n _i 

The right hand side of F^H*) is positive, and decreases as i increases, 
therefore F,(H*) needs to be positive and decreasing which is satisfied 
by 


(2.4) 
and 

(2.5) 


Qn ^ lb 


ri > r 2 > ■ ■ ■ > r r 


Inequalities (12. 4 p and (12.5p can be combined to give 
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Proposition 2.1. E* exists if and only if \2.2\) . | ITfj) and \2.5\) hold. 
In fact, 


(2.7) 


HI = 


H* = - 


ei 




3 > 1 , 


Pf = rj - r j+ 1 , j < n, P* = r n — Q r 


Furthermore, the above also implies the existence of a unique equilib¬ 
rium E i with all components positive except for P n = 0. In fact, 


(2.8) Hi = H* n + H) = H*, 1 <j<n, 

a n J 

P] = P*, 3 < n, Pi = 0. 

Remark 2.2. (12.61) implies the existence of an unique family of equilib¬ 
ria El with Hj , Pj = 0, j > k described by (12.71) . but with Q k replacing 
Q n . Another family of equilibria, E' k , exists with Hj = 0, j > k and 
Pj = 0, j > k described by (j2.8[) . but with Q k replacing Q n . There are 
many other equilibria but we have no need to enumerate all of them. 


Hereafter, we assume without further comment that (12.4p . and (12.51) 
hold. 

We write H i oo = lim inf H t (t) and H°° with limit superior in place 

t—too 

of limit inferior. 

Proposition 2.3. (a) If (^ Hj(t))°° < j- then Pi(t) — > 0. 

j<i 

(b) Ifi < j, Pi( 0) > 0, and if (H i+1 + H i+2 H-f Hj)°° < j- - P 

then Pj(t ) —>■ 0. 

(c) Ifi < j , Hf 0) > 0, and if ( P t + P i+1 H-f Pj-i)°° < ( r i ~ r j) 

then Hj(t ) —> 0. 

Proof, of (a): The equation for P, implies that 


f t \o g p--‘ 


j<i 


1 




n 

If (X) Hj(t))°° < j- is false, then Pi —» oo, a contradiction to Pi being 

i<i 

bounded. Assertion (a) is transparent. 
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Proof of (b): If i < j, P,(0) > 0, and if (H l+ \ + H i+2 + • • • + Hj)°° < 
-— — then 

Cj €-i 


1 _ 

n,A e,- 


d Pf iH 

P og 


1_ 

riA e.A 


p '“J °3 

3 


= _P[ _ 

n j e j P j 

= - 1 - {H i+ 1 + Hi _|_2 + • • • + Hj) > e, f > T 

e* Cj 


JID 

for some e, T > 0. Therefore, ‘ i —> oo, and since Pi, Pj are bounded, 

p n j e 3 

Pj(t ) o. 

Proof of (c): assume that Hf 0) > 0 and (Pj + P* + i H— • + Pj_i)°° < 
(■ r i — Tj), then 


± [ae ML H '< H '= 

dt 8 Hj(t) Hi Hj 

— { r i ~ r j ) — (Pi(t) + Pi+l{t) + ‘ ‘ ‘ + Pj-i(t)) 


It follows that ^ —> oo, and since HHj are bounded, Hj{t) —> 0. □ 

If there are no virus present, then host H\ drives the other hosts to 
extinction. 

Lemma 2.4. If Pi = 0,1 < i < n, Pi(0) > 0 then Hi —» 

Proof. Since P, = 0, //, +1 —* 0 by Proposition 12.31 (c) for 1 < i < n. 

n 

Therefore Ve > 0, 3T > 0 such that Vt > T, a jHj(t) < e. Then for 

3 =2 

t > T,H[ > H\ (r'! — aiHi — 2e). Therefore Pi i0O > ri ~ 2<: and since 
e > 0 is arbitrary, H i iOQ > —. On the other hand, H[ < H\{ri — a\H\), 
so Hr < a. Therefore PW ^ □ 

Now we show that Hi persists if initially present regardless of who 
else is around; similarly, Hi and V\ persist if initially present regardless 
of which other host and virus are present. 

Proposition 2.5. (a) If H\ (0) > 0, then H^° > T. 

(b) If H\ (0) > 0 and Pi(0) > 0 , then 

Hr > —, P?° > minln - r 2 , riei ~ ai }. 

ei ei 
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Proof, of (a). Assume the conclusion is false. Then P\ —y 0 by Propo¬ 
sition 12.31 (a). 

If Pj(0) = 0 for all i, then Hi(t) —> > A by Lemma [2~P and (12.6p . 
so we suppose that P;(0) > 0 for some i. Let k denote the smallest 
such integer i for which Pj(0) > 0. 

If k — 1 , then, as noted above, Pi —> 0 and so H 2 — » 0 by 
Proposition 12.31 (c). Then P 2 —> 0 by Proposition 12.31 (a) or (b). 

If k — 2, then P\ = 0 so H 2 —> 0 by Lemma [2.31 (c) since Pi and P 2 
share the same virus. Since (Pi + P 2 )°° = Pf° < A < A, it follows 
that P 2 —> 0 by Proposition 12.31 (a). Now we can use Proposition ^. 31 (c) 
to show H 3 —> 0 and then Proposition 12.31 (a) or (b) to show P 3 —>■ 0. 

If k > 2, then P x = P 2 = ■ ■ ■ = Pk-i = 0 and Pfc(0) > 0. As 
Pi, ■ • ■ , P/,._] share the same virus, then //, = 0 or P ?; —» 0 for 1 < i < 
k — 1 by Proposition 12.31 (c). H k —» 0 by Proposition 12.31 (c). Then, 
P fc —> 0 by Proposition 12.31 (a) since (^ Pj)°° = < A < A_ g Q 


Pfc + i —> 0 by Proposition 12.31 (c). Proposition 12.31 (a) or (b) implies 
that Pfc+i 0. 

We see that for all values of k, P 2 , • • • , H k + 1 —> 0 and Pi, • • • , P k +\ —>■ 
0. Successive additional applications of Proposition 12.31 (a) or (b) and 
(c) then imply that P 2 , • • • , H n —» 0 and Pi, • • ■ , P n —* 0. But, then 
for all e > 0, there exists T > 0 such that 


P(/(Pi) > ri - e - aiPi, t > T. 


This implies that Pf° > ^ > A, by (12.41) . a contradiction. This com¬ 
pletes the proof of the first assertion. 

Proof of (b): Now, suppose that Pi(0) > 0,Pi(0) > 0 and Pf° < 

ri - r 2 . 

Proposition 12.31 (c) implies that H 2 —> 0. By Proposition 12.31 (b), 
P 2 —> 0. Applying Proposition 12.31 (c) with i — 1 and j = 3, as 
(Pi + P 2 )°° = P{° < ri — r 2 < ri — r 3 , we conclude that —> 0. Then, 
Proposition 12.31 (b) implies that P 3 —y 0. Clearly, we can continue 
sequential application of Proposition 12.31 (b) and (c) to conclude that 
Hi, Pi —> 0 for i > 1. Then we use that 




a l 

e l n \ 


Hi Pieini 

di „ . 

ri-Pi — terms that go to zero 



to conclude that P?° > riei ai . 

1 — ei 


□ 
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Lemma 2.6 (Lemma 1.2 [6]). Let x{t) be a bounded positive solution 
of the Lotka- Volterra system 

n 

x\ = Xilri + a,ijXj), 1 < i < n 

3 = 1 

and suppose there exists k < n and m, M , 5 > 0 such that m < Xi(t) < 

M, 1 < % < k, t > 0, Xk+i(t) <5, t > 0, and Xj(t ) —* 0 for j > k + 1. 

Suppose also that the k x k subsystem obtained by setting x 3 — 0, j > k 
has a unique positive equilibrium p = (pi,p 2 , • • • ,p k )- Then 

i r T 

liminf — / xAt)dt — pi + 0(8), 1 < i < k. 

T^oo T J 0 

The same expression holds for the limit superior. 

Theorem 2.7. Let 1 < k < n. 

(a) There exists e k > 0 such that if Hi (0) > 0, 1 < i < k and 
Pj( 0) >0, 1 < j < k — l, then 

Hi ,oo > e k , 1 < i < k and P jj00 > e k , 1 < j < k - 1. 

(b) There exists e k > 0 such that if Hf 0) > 0, Pj(0) >0, 1 < i < k, 

then 

Hi t oo > e k , P i,oo > e k , 1 <i < k. 

Proof. We use the notation [Hf\ t = j f* Hi(s)ds. Our proof is by math¬ 
ematical induction using the ordering of the 2 n cases as follows 

(a, 1) < ( b, 1) < (a, 2) < (6, 2) < • • • < (a, n) < ( b , n) 

where (a, k) denotes case (a) with index k. 

The cases (a, 1) and ( b , 1) follow immediately from Proposition 12.51 
and Corollary 4.8 in [4j with persistence function p = min{//i, P \} in 
case (6,1). The latter result says that weak (lirnsup) uniform persis¬ 
tence implies strong (liminf) uniform persistence when the dynamical 
system is dissipative. 

For the induction step, assuming that (a, k) holds, we prove that 
(6, k ) holds and assuming that (6, k ) holds, we prove that (a, k + 1) 

holds. 

We begin by assuming that (a, k) holds and prove that (6, k) holds. 
We consider solutions satisfying Hi( 0) > 0,Pj(0) > 0 for 1 < i < k. 
Note that other components Hj( 0) or Pj( 0) for j > k may be positive 
or zero, we make no assumptions. As (a, k ) holds, there exists e k > 0 
such that Pj,oo > efc, 1 < i < k and Pj i00 > e k , 1 < i < k — 1. We need 
only show the existence of 5 > 0 such that P k)00 > S for every solution 
with initial values as described above. In fact, by Corollary 4.8 in [4], 
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weak uniform persistence implies strong uniform persistence, it suffices 
to show that Pg° > 6. 

If P™ < r k — r k+ i, then H k+i — * 0 by Proposition 12.31 (c). Then, 
by Proposition 12.31 (b), P k +\ —> 0. Clearly, we may sequentially apply 
Proposition 12.31 (b) and (c) to show that H } —* 0, Pj —> 0 for j > k + 1. 

If there is no 6 > 0 such that P£° > 6 for every solution with initial 
data as described above, then for every S > 0, we may find a solution 
with such initial data such that P^° < 5. By a translation of time, 
we may assume that P k (t) < S, t > 0 for 0 < 5 < r k — r fc+1 to be 
determined later. Then Hj , Pj —> 0, j > k + 1. Now, as (a, k) holds, we 
may apply Lemma 12.(11 The subsystem with II, = 0. k+1 < i < n and 
Pi = 0, k < i < n has a unique positive equilibrium by Proposition ^. 11 
See Remark 12.21 The equation 


p> 

r k 

f A' I k 11 k 




1 


implies that 


b°g 


^ - £1^]. 

p e k n k 


(0) j<k 

By (12.8p and Lemma [2.61 we have for large t 


1 


Ei^i< - 

j<k 


1 

e-k 


E*J 


-t + 0 ( 6 ) 

&k 


— + q-- + 0(5) 

^k —1 &k 


where q = —(r k —Q k - 1 ) + —-— > 0. On choosing 5 small enough and 

a k &k — 1 &k 

an appropriate solution, then [Hj]t —- > q/2 for large t, implying 

j<k ek 

that P k -P +cxo, a contradiction. We have proved that (a, k ) implies 
(6, k). 

Now, we assume that (6, k) holds and prove that (a, k + 1) holds. We 
consider solutions satisfying 

Hi{ 0) > 0, Pj(0) > 0 for 1 < i < k and H k+ 1(0) > 0. As (6, k) holds by 
assumption, and following the same arguments as in the previous case, 
we only need to show that there exists 5 > 0 such that H^_ 1 > 5 for 
all solutions with initial data as just described. 

If H^_ 1 < then P k +i —> 0 by Proposition 12.31 (b) and then 

H k+ 2 —> 0 by Proposition 12.31 (c). This reasoning may be iterated to 
yield H, -P 0, k + 2 < i < n and P % -P 0, k + 1 < i < n. 

If there is no 5 > 0 such that H^_ 1 > 6 for every solution with 
initial data as described above, then for every 5 > 0, we may find a 
solution with such initial data such that H^_ 1 < 6. By a translation of 
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time, we may assume that H k +\(t) <5, t > 0 for 0 < 5 < — — 

to be determined later. Then Hj,Pj — * 0, j > k + 2 and Pk+i —> 0. 
Now, using that (6, k) holds, we apply Lemma 12.61 The subsystem 
with Hi — 0, Pi — 0 k+1 <i <n has a unique positive equilibrium by 
Proposition 12.11 See Remark 12.21 The equation for H k+1 is 


H' k+ 1 


Hk+i 

Integrating, we have 


= r k+ 1 


a jHj a jHj 

j=l j=k +1 


n 


E a 

j=fc+i 


ilog*±lW 


fc+1 


j=i 


By Remark 12.21 and Lemma 12.61 we have that for all large t 
By (12.7p and Lemma 12.61 we have that for all large t 


k k 

Y a j[Hj] = Y a i H ] + °( 5 ) = Qn + 0(5) 

3 = 1 3 =1 

since H k+ i(t) < 5, [H k +i]t — 0(h). Now, Q n > 0 so by choosing 5 
sufficiently small and an appropriate solution, we can ensure that the 
right hand side is bounded below by a positive constant for all large 
t, implying that H k+ i(t) is unbounded. This contradiction completes 
our proof that (6, k) implies (a, k + 1). Thus, our proof is complete by 
mathematical induction. □ 

Corollary 2.8. For every solution of (12.ip starting with all compo¬ 
nents positive, we have that 

(2.9) j [ Hi(s)ds -A H *, \ [ Pi(s)ds^P* 

t Jo t Jo 

where H*,P* are as in \2.7L 

For every solution of (12.ip starting with all components positive ex¬ 
cept 

P n (0) = 0, we have that 

(2.10) - f Hi(s)ds -A Hi - f P,{s)ds P} 

t Jo t Jo 

where H\,P} are as in (12. 8p . 


Proof. This follows from the previous theorem together with Theorem 
5.2.3 in 0. □ 
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3. Global Stability for NIN in a special case 
Positive equilibrium E* exists so we can write the system as 


(3.1) 




0=1 


P' = e,n. 


Pi (l>i - h;) j , 


j>i 

1 < i < n. 


Let U(x,x*) = x—x* — x* logx/x*, x,x* > 0, be the familiar Volterra 
function and let 

v=Y, c ' u[H ^ + E d - u{p ^ p t) 


where ci, ■ • • , c n and di, ■ ■ ■ ,d n are to be determined. 

Then the derivative of V along solutions of (13.11) . V, is given by 


V 


Y.dieMPi- pnY.Wi - H i) 

i j<i 


i j>i 


We aim to choose parameters so that the last two terms cancel each 
other. The second summation may be rewritten as ^(Pj— P*) Y2 c j{Hj — 

i j<i 

H*) so that the last two sums may be combined as ^(P—P*) — 

i j<i 

Cj)(Hj — H*). It vanishes if Vi, djejUj — Cj — 0, j < i. Taking i — n, 
we see that the Cj must be identical so Cj = a for all j for some a > 0 
and di = a/e^ru. Therefore, in this case, we have 

V = - a(H, - H*)j a t (H, - if*)J 


If, in addition, a.j = a for all j, then we have 



Theorem 3.1. Assume that holds arid aj = a > 0, 1 < j < n 
in (13.11) . Then E* is globally asymptotically stable relative to the open 
positive orthant ofM. 2 ^. 
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With the same assumptions, but Hi, Vj replacing H*, V* and addition¬ 
ally P n = 0, El is globally asymptotically stable relative to the open 
positive orthant 

Proof. We first note that since V(H(t),P(t )) < 14(0)), t > 0 

for every positive solution of (13. ip . each component is bounded above 
and below: 0 < p < x(t ) < P, t > 0, where x = H{,Pj and p, P may 
depend on the solution. 

Consider a positive solution of (13.ip . By LaSalle’s invariance princi¬ 
ple, every point in its (invariant) limit set L must satisfy 22 Hi = 22 H* 

i i 

since L C {(H, V) : V = 0}. Since V{x) < V{H( 0), P(0)) for all x e L, 
L belongs to the interior of the positive orthant and it is bounded away 
(but maybe not uniformly) from the boundary of the orthant. We now 
consider a trajectory belonging to L; until further notice, all consider¬ 
ations involve this solution. Notice that this solution satisfies 

(3.3) H[ = H, ng(P; - P,)J 

(3.4) p; = ei n,p, - H’)j 

From (13.4p . we see that P' n = 0 so P n {t) is constant. Then, H' n = 
H n (P* — P n ) so H n (t) is either converging exponentially fast to zero, 
blowing up to infinity, or identically constant depending on the value 
of P n . The only alternative that is consistent with L being invariant, 
bounded, and bounded away from the boundary of the orthant is that 
H n (t) is constant and that P n = P*. As we use a similar argument 
repeatedly below, we refer to it as our standard argument. 

Since H n is constant and JT Hj is constant, equal to JT H*, then so 
is 22i < n -1 Hi a constant. But now we face the same dilemma as above 
with the equation (I3.4p with i = n — 1 since the sum in parentheses 
is constant. By our standard argument, the only alternative is that 
this constant is zero, i.e., that Hi = X^<n-i-^T anc ^ Pn-iif) is 

constant. The former implies that 

= E = H 'r 

i i<n —1 

Suppose that 1 < k < n and that H,ft) = H*, Pjft) = P*, i > k, 
hold. We claim that = Hl_ v P k _i(t) = P k _ v As P k (t) is 

constant, f!3.4p implies that -< fc Hj(t) = J2j<k Hj anc ^ since H k = H£, 
it follows that Ylj<k-iHj(^) = Ylj<k-iHj- N°fi ce that if k — 2, then 
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the latter gives that H\ = H{. Now from (13.41) . P/,_ 1 (t) = 0 so Pk-i{t ) 
is constant. This implies, by (13.3p and our standard argument, that 
P^._ 1 = 0 and Pk-\(t) = Pf_ v If k — 2, we are done: Hi = H*, Pi = 
P*- If k > 2, then J2j<k-2 H A t ) = ^2j<k -1 H A t ) ~ H k-iif) is constant 
so from (13.41) and our standard argument we conclude that P[_ 2 = 0 
and that Hj(t) = )T) ;<fc _ 2 H*. The latter implies that 

ft-. = v Hj- Y, H i= E H i~ E H i= 

j<k-1 j<k -2 i<fc-2 

This completes our proof of the claim. By induction, we conclude that 
Hi(t) = H*, Pi(t) = P*, 1 < i < n, i.e., our solution is identical to 
E*. Since we considered an arbitrary solution starting at a point of 
L , it follows that L = {P*}. As our chosen solution was an arbitrary 
positive solution, we have established the result. 

The arguments are nearly identical for the case. From (13.31) . H' n = 0 
since P n = 0, therefore the standard argument starts at n — 1 instead. 

□ 

4. One-to-One Infection Network 

M — I in the one-to-one infection network so the equations then 
becomes: 


(4.1) 


H[ = II, ^ ~ HiPi 

P'i = eiTiiPi K Hi - A ] , 1 < i < n. 


The principle equilibrium for the one-to-one infection network are 
now described. 


Proposition 4.1. There exists an equilibrium E* with Hi and Pi pos¬ 
itive for all i if and only if the following inequality holds: 


(4.2) Q n < r v 1 <j<n, Q n = ^ —• 

i=l 

In fact, 

(4.3) H * = -, j > 1, 

e 3 

Pj = Tj ~ Q n , j > 1 . 

The positive equilibrium E* is unique. 
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We also note the existence of a unique equilibrium E\ with all com¬ 
ponents positive except for P n — 0 , given by 

H) = 

(4.4) Hi = K + H, 

&n 

Pj = Pj ~ Pn = r 3 ~ j < n, 

provided that r n < rj,j ^ n . 

Remark 4.2. We will assume hereafter that 

(4.5) r L > r 2 > ■ ■ ■ > r n . 

This hypothesis ensures the existence of a family of equilibria El , E' k , 1 < 
k < n, characterized as follows. El with Hj , Pj = 0, j > k is described 
by (14.31) but with Qk replacing Q n . E\ satisfies Hj = 0 , j > k and 
Pj = 0, j > k described by (14.41) but with Qu replacing Q n . 

Proposition 4.3. (a) If H.ft) 00 < A then Pi(t) —> 0. 

(b) If i < j, Hi( 0) > 0, and (Pi — Pj)°° < r* — rj, then Hj(t ) —>■ 0. 

Proof, of (a). The equation for Pj implies that 


4 log />■■•■ =*,(*)-4 


If Hift) 00 < T then Pj -)• 0. 


Proof of (b). Assume that i < j, Pj(0), Hj( 0) > 0 and (Pj — Pj 


< 


rj. As 


d . Hi(t) 
dt X ° g Hj(t) 


H[_El 

Hi H, 


= (n-Tj) - (Pi(t) -Pj(t)), 


it follows that ^ —> oo, which by the boundedness of Hi,Hj, implies 
that Hj(t ) —y 0. □ 


Lemma 4.4. If Pi = 0,Pi(0) > 0 then Hi —> 

Proof. Since Pi = 0, H t —> 0 by Proposition 12.31 (b) for 1 < i < n. 

n 

Therefore, Ve > 0, 3 T > 0 such that Vf > T, a jHj(t) < e. Then for 

3 = 2 

t > T, H\ > H\ (rj — aiHi — 2e). Therefore Pi i00 > ri a ) 2e and since 
e > 0 is arbitrary, Pi )00 > —. On the other hand, H[ < Hi(ri — aiHi), 
so H?° < a. Therefore □ 

Proposition 4.5. (a) If H\ (0) > 0, then Hf° > A. 
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(b) If H i(0) > 0 and -Pi(O) > 0, then 


1 

ei 


Hf° > —, Pf° > min{ri — r 2 , 


rie i — a i 
ei 


}• 


Proof, of (a): Assume the conclusion is false. Then Pi —> 0 by Propo¬ 
sition 14.31 (a). Then H 2 —:> 0 by Proposition 14.31 (b). Therefore by 
sequential applications of Proposition EO (a) and (b), we can conclude 
that Hi , Pi —> 0, for i > 1. But, then 

H[l(HP > ri — e — aiHi > r n — e — aiHi > — — e — aiHi, t >T 

ei 

for some e > 0 and T > 0 (recall that r t > r n from (14.5p and ^ < 
7T from (14.2p ). This implies that Hf° > A., a contradiction. This 
completes the proof of the first assertion. 

Proof of (b): Now, suppose that Pi(0) > 0, P\ (0) > 0 and Pf° < 
7T — r 2 . Then (P\ — P 2 )°° < Pf° < T\ — r 2 , therefore Proposition 14.31 (b) 
implies that H 2 —> 0. Then P 2 —y 0 by Proposition 14.31 fa). Then (P 2 — 
-P 3) 00 < Pf° < r 2—^3 therefore Proposition 14.31 (b) implies that H% —> 0. 
Clearly, we can continue sequential applications of Proposition 14.31 (a) 
and (b) to conclude that Hi, P, : —* 0 for i > 1. 


| log HiP™ 


H[ aiPj 

Hi + Pieini 
ai 

ri -Pi — terms that go to zero 

ei 


to conclude that P?° > riei ai . 

1 — ei 


□ 


Theorem 4.6. Let 1 < k < n. 

(a) There exists > 0 such that if Hi( 0) > 0, 1 < i < k and 
Pj( 0) >0, 1 < j < k — 1, then 

Hi ,oo >e k , 1 < i < k and P j)00 > e k , 1 < j < k - 1. 

(b) There exists e k > 0 such that if Hf 0) > 0, Pj(0) >0, 1 < i < k, 
then 

Hi ,00 > ejfc, P, 1,00 >e k , 1 < i < k. 


Proof. We use the notation [Hft = j f ( ' Hi{s)ds. Our proof is by math¬ 
ematical induction using the ordering of the 2 n cases as follows 

(a, 1) < (6,1) < (a, 2) < (b, 2) < • • • < (a, n ) < (6, n) 

where (a, k) denotes case (a) with index k. 

The cases (a, 1) and ( b , 1) follow immediately from Proposition 14.51 
and Corollary 4.8 in [4J with persistence function p = miri{/p, P\ } in 
case (6,1). 
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For the induction step, assuming that (a, k) holds, we prove that 
(6, k) holds and assuming that (6, k) holds, we prove that (a, k + 1) 
holds. 

We begin by assuming that (a, k ) holds and prove that (b, k) holds. 
We consider solutions satisfying H^O) > 0, P,;(0) > 0 for 1 < i < k. 
Note that other components Hj(0) or Pj( 0) for j > k may be positive 
or zero, we make no assumptions. As (a, k) holds, there exists e k > 0 
such that H ioo > e k , 1 < i < k and P i)OQ > e k , 1 < i < k — 1. We 
need only show the existence of 8 > 0 such that P k:00 > 8 for every 
solution with initial values as described above. In fact, by the above- 
mentioned result that weak uniform persistence implies strong uniform 
persistence, it suffices to show that P£° > 5. 

If Pg° < r k — r k+ 1 , then H k +\ —» 0 by Proposition 14.31 (b). Then, by 
Proposition 14.31 (a). P k+ \ —> 0. Clearly, we may sequentially apply 
Proposition 14.31 fbf and (a) to show that Hj —> 0, Pj —)■ 0 for j > k + 1. 

If there is no 5 > 0 such that P£° > 8 for every solution with initial 
data as described above, then for every 5 > 0, we may find a solution 
with initial data such that P£° < 8. By a translation of time, we may 
assume that P k (t) < 8, t > 0 for 0 < 8 < r k — r fc+1 to be determined 
later. Then Hj,Pj —> 0,j > k + 1. Now, as ( a,k ) holds, we may 
apply Lemma 12.61 The subsystem with Hi — 0, k + 1 < i < n and 
Pi — 0, k < i < n has a unique positive equilibrium by Proposition l4.ll 
See Remark 14.21 The equation 



Pk&kN'k Gfc 


implies that 


By (14.4p and Lemma [2.61 we have for large t 

[H k \ t - = Hi -b 0(5) = Hl~\ ——-b 0(5) = — + 0(8) > 0 

ek &k C^k 

Implying that P k — > +cxo, a contradiction. We have proved that (a, k ) 
implies (b,k). 

Now, we assume that (6, k ) holds and prove that (a, k + 1) holds. We 
consider solutions satisfying 

Hi{ 0) > 0, Pi{ 0) > 0 for 1 < i < k and Lffc + i(0) > 0. As (6, k) holds by 
assumption, and following the same arguments as in the previous case, 
we only need to show that there exists 8 > 0 such that H^_ 1 > 8 for 
all solutions with initial data as just described. 




1 

D e k n k 

r k 


( 0 ) 
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if m 


< 


-, then P k _|_i —y 0 by Proposition 14.31 (a) and then 


H k+ 2 —> 0 by Proposition 14.31 (b). This reasoning may be iterated to 
yield Hi —> 0, k + 2 < i < n and Pi 0, k + 1 < i < n. 

If there is no 8 > 0 such that H ^ > 8 for every solution with 
initial data as described above, then for every 8 > 0, we may find a 
solution with such initial data such that H^f 1 < 8. By a translation 
of time, we may assume that H k+ i(t) <8, t > 0 for 0 < 8 < —— 
to be determined later. Then Hj,Pj —> 0 ,j > k + 2 and P k+ i —> 0. 
Now, using that (6, k ) holds, we apply Lemma 12.61 The subsystem 
with Hi = 0, Pi = 0 k + 1 < i < n has a unique positive equilibrium by 
Proposition 14.11 See Remark 14.21 The equation for H k+ \ is 


H' k+ 1 


H, 


k +1 


r k+i ^ a jiH] 'y a jHj p k +\ 
j =1 j=k +1 


Integrating, we have 


1 g t+ i( i> 
t °H t+ ,(0) 


fc +1 

= y ^ a j[Hj\t + o(i/t) 

j= 1 


By (14.3p and Lemma 12.61 we have that for all large t 

k k 

y ^ cij[Hj] = y ^ cijH* + 0(8) = Q n + 0(8). 
i = 1 i =1 

Since H k+ i(t) < 8, [H k+ i) t = 0(<5). Now, Q n > 0 so by choosing 8 
sufficiently small and an appropriate solution, we can ensure that the 
right hand side is bounded below by a positive constant for all large 
t , implying that H k+ i(t) is unbounded. This contradiction completes 
our proof that (6, k) implies (a, k + 1). Thus, our proof is complete by 
mathematical induction. □ 


Corollary 4.7. For every solution of (14.ip starting with all compo¬ 
nents positive, we have that 

(4.6) y [ Hi(s)ds —>• H *, i [ P t (s)ds -+ P* 

t Jo t Jo 

where H*,P* are as in (14.31) . 

For every solution of (14.11) starting with all components positive ex¬ 
cept P n (0) = 0, we have that 

(4.7) - f Hi(s)ds —>• Hj, - f Pi(s)ds —>■ P} 

t Jo t Jo 

where Hj,P- are as in (14.4p . 
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Proof. This follows from the previous theorem together with Theorem 
5.2.3 in 0. □ 

5. Global dynamics for the one-to-one network 
Using the positive equilibrium E*, we can write the system as 


(5.1) Hi = H, -h s ) + p;-P. 

P[ = eiUiPi (Hi - H*) , 1 < i < n. 

As before, we define 

v = Ci U(Hi,H*) + J2 diU(P i: P*) 

i i 

where ci, ■ ■ ■ ,c n and d\, ■ ■ ■ ,d n are to be determined. 

Then the derivative of V along solutions of (15.11) . V, is given by 

v = - (£- H’fj aj (H s - H*)J - £ Ci (H, - H*)(Pi - P-) 
+ J2 d <e,n,(P l - H‘) 

i 

Letting c t = a,i and d* = - a - causes the last two summations to cancel 
each other out. Therefore in this case we have 



Below, we use the notation (H(t), P(t)) for the 2n-vector solution 
,if n (t),Pi(t),..- ,P n (t)). 

Theorem 5.1. The co-limit set of a positive solution of ( 15 . 11 ) is either 
E* or it consists of non-constant entire orbits, (H(t), P(t)), satisfying 
all of the following: 

(a) E”=i Oi W) = E"., Oi h;, t 6 E. 

(b) YYUm^ is independent of t. 

(c) Vi, (Hi(t), Pi(t)) is a positive solution of the conservative planar 
system 

= H^Pf-Pf) 

= e l n l P t (H t - H*). 


(5.3) 


P! 
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Similarly, the co-limit set of a solution of (15.1ft with all components 
positive except P n = 0 is either E t or it consists of non-constant en¬ 
tire orbits, as in the previous case, on the hyperplane ^” =1 a^H^ft) = 
Y^i=i a iH} with H n (t ) = Hf and with n”=i' E i (t) ai ^ eini independent of 
t. Furthermore, for 1 < i < n, (Hft), Pi(t)) satisfies (15.3p but with 
Hi , V} replacing H*, V*. 

Proof. We first note that since V(H(t), P(t)) < V(H(0),V(0)),t > 0 
for every positive solution of (15. ip . each component is bounded above 
and below: 0 < p < x(t ) < P,t > 0, where x = H i: Pj and p, P may 
depend on the solution. 

Consider a positive solution of (15.1ft . By LaSalle’s invariance prin¬ 
ciple, every point in its (invariant) limit set L must satisfy JA ajPj = 
JA a-iH* since L C {( H , V) : V — 0}. As in the NIN case, L belongs 
to the interior of the positive orthant and it is bounded away from 
the boundary of the orthant. We now consider a trajectory belonging 
to L ; until further notice, all considerations involve this solution. No¬ 
tice that this solution satisfies (15.3p . Thus on L, the system decouples 
into n independent planar conservative systems, the positive solution of 
which is either periodic or is the positive equilibrium. See e.g. section 
2.3 of [2|. Notice that fP ll “'ff, = — #*) = 0, consequently 

K=,m ai/eini i s independent of t. 

If E* G cc-limit set, then E* = cu-limit set, since E* is stable. Conse¬ 
quently, if E* fz L, then at least one of the (H,, Pf) must be a non-trivial 
periodic orbit. 

The arguments are nearly identical for the case that the solution 
satisfies P n = 0 and other coordinates positive. Liapunov function V 
differs from the previous one only in that the sum goes from one to 
7i—l in the second summation and Hj, V} replace H*, V*\ the choice 
of the Ci and di are as before. (15.2p is changed only in that superscript 
f replaces *. 

We only note that the counterpart to (15.311 for i = n reads H' n = 0. 
As JA ai(Hi — Hj) = 0 on the limit set and since any positive periodic 
limiting solution must satisfy J Q T Hidt = Hj, it follows that H n = Hf. 

□ 


In the special case that n = 2, since Hi (Pi) can be expressed in 
terms of H 2 (P 2 ), on {(P, V) : V = 0}, every solution in L is periodic 
(possible constant). 
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